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Abstract 

Absorbing boundary condition approach to nuclear breakup reactions is investigated. A key 
ingredient of the method is an absorbing potential outside the physical area, which simulates 
the outgoing boundary condition for scattered waves. After discretizing the radial variables, the 
problem results in a linear algebraic equation with a sparse coefficient matrix, to which efficient 
iterative methods can be applicable. No virtual state such as discretized continuum channel needs 
to be introduced in the method. Basic aspects of the method are discussed by considering a nuclear 
two-body scattering problem described with an optical potential. We then apply the method to the 
breakup reactions of deuterons described in a three-body direct reaction model. Results employing 
the absorbing boundary condition are found to accurately coincide with those of the existing 
method which utilizes discretized continuum channels. 
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I. INTRODUCTION 



In nuclear reactions with weakly-bound projectiles, breakup processes are significant not 
only for their prominent cross sections as a final reaction product but also for their influ- 
ences on other reaction channels such as elastic scattering. The breakup processes have 
been investigated originally for the reactions involving weakly-bound stable nuclei, such as 
deuterons 0] and ^'''Li 0. In the last decade, interests on the breakup reactions have 
explosively increased because of the discovery of the halo nuclei around the neutron drip-line 
and the studies of their reactions 

A variety of theories have been developed to describe the breakup processes . In this ar- 
ticle, we consider reactions in which the weakly-bound projectile is described as a composite 
two-body system and the whole reaction is described as a three-body problem of the projec- 
tile and the target. At sufficiently high incident energies, one may assume that the projectile 
internal motion is much slower than the projectile-target relative motion. Then, the projec- 
tile internal motion may be frozen during the reaction. This is the adiabatic approximation, 
which has been developed originally for the deuteron reactions [^]. One may also adopt 
the eikonal approximation for the fast projectile-target relative motion. A combined use of 
the eikonal and adiabatic approximations has provided especially useful descriptions for the 
reactions of halo nuclei 0. 

For a more accurate description, a fully quantum mechanical treatment of the breakup 
processes has been developed in a coupled channel framework. This is known as the coupled 
discretized continuum channels (CDCC) method ^, In this method, the continuum 
excited states of the projectile are incorporated by discretizing the momentum variables. 
The CDCC method has been successful in describing the reactions of deuterons H and 
^'''Li 0, and has been recently applied to the reactions of halo nuclei 0. 

In the present article, we discuss an alternative, fully quantum mechanical method: the 
method of absorbing boundary condition (ABC). This method was developed in the field of 
chemical reactions ^, [l^. The key ingredient is an introduction of the absorbing potential 
outside the physically relevant area. This absorbing potential allows us to treat scattering 
problems with the vanishing boundary condition. We show that the ABC method provides 
a useful description for the breakup reactions in which the final scattering channels include 
three-body continuum states. 
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The present article is organized as follows: In Section II, we illustrate the ABC method 
in the simplest example, a two-body scattering problem with a nuclear optical potential. In 
Section III, we formulate the ABC method for three-body scattering problems. In Section 

IV, we present results of the ABC method applied to the breakup reaction of rf+^^Ni at 
the incident deuteron energy of 80 MeV. The results are then compared with those by the 
CDCC method. The summary is given in Section V. 



II. ABSORBING BOUNDARY CONDITION METHOD: A POTENTIAL SCAT- 
TERING 

We first explain the ABC method considering the simplest example: the scattering of a 
particle of mass m from a potential V{y). We separate the potential V{v) into two parts, 

y(r)=yo(r) + AT/(r). (1) 

The potential Vo(r) is spherical and may be of long-range . We assume that it is easy to 
construct the scattering solution ■0^'^^(r) for the potential Vq{v). The potential Ay(r) may 
not be spherical but must be of short-range , Ay(r) = for the spatial region beyond a 
certain radius rc, |r| > Tc- 

The scattering wave function can be written in a form 

^(+) (r) = (r) + AV'^^) (r) . (2) 

The function AiP'^^\y) is expressed by 

A^W(r) = j dr'G^^\vy )/\V{Y')ij^+\v') = J dr'G'W(r, r')Ay(r')4^\r'), (3) 

where Gg^^r, r') and G*^"*") (r, r') are the Green's functions with the potentials Vo{r) and 
V{r), respectively, 

= , G W = . (4) 

° E + ie-n-Vo E + ie-T^-V ^ ^ 

Here, is the kinetic energy operator. A positive infinitesimal number e specifies the 

outgoing-wave boundary condition. For the special case of K) = with AV = V, 4^V) 

becomes a plane wave and we have an analytic expression for G'q~\ 

Now, let us explain the ABC method to construct the wave function A'0*^+)(r). The basic 

trick is that the positive infinitesimal number e is replaced by a finite, space-dependent 



function e(r). In order to simulate the outgoing boundary condition, we introduce the 
function e(r) which is positive and finite in the spatial region, Vc < r < rc + Ar. We impose 
the vanishing boundary condition at r = + Ar for Aip^~^\v). If the absorbing potential 
works ideally, only the outgoing waves are allowed in the spatial region around r = Vc- If 
we regard e(r) as a part of the Hamiltonian, namely Ti. = Tr + V — ie{r), one may regard 
this replacement as an addition of the absorbing potential, —ie{r), to the Hamiltonian, 
H = Tr + V. 

Employing the absorbing potential e(r), we can rewrite Eq. @ for A'ip^^\r) as the fol- 
lowing linear inhomogeneous equation, 

{E + ie{r) - Tr - y)A^W(r) = Ar(r)4^V)> (5) 

where, as mentioned above, the scattered wave Atp^~^\r) should satisfy the vanishing bound- 
ary condition at r = rc + Ar. The right hand side of this equation is spatially localized since 
we assume that the potential AV{r) vanishes in the spatial region r > Tc- 

The ABC approach thus allows us to obtain the scattered wave function in the spatial 
region of r < by imposing the vanishing boundary condition at r = + Ar. We should 
note that no asymptotic form is required to solve the scattering problem in this procedure. 
The elastic scattering amplitude /(k) can be obtained from the wave function in the spatial 
region where the potential AV{r) exists, 

/(k) = /o(k) - ^ / rfr4k^*(r)A\/(r)^W(r), (6) 

where /o(k) is the scattering amplitude with the potential Vo(r) and 'ipoiii''^) incoming- 



wave solution in the potential Vo(r) with the incident momentum k (/c = a/ 2mE / h and 
k = k/A;). Here, the wave functions are normalized as |'?/'(r)| — 1 at r — oo. 

The accuracy of the calculated solution depends on the quality of the absorbing potential. 
One should choose the function e(r) carefully to make the reflected waves as small as possible. 
In this respect, the absorbing potential with linear dependence on the coordinate has been 
tested extensively [|1^, 0. We parameterize the potential as 

(r < Tc) 

te(r] = < , (7) 

-iW^hs'^ (r^ < r < rc + Ar) 
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where Wabs, and Ar are positive constants, representing the strength, radius, and thick- 
ness of the absorbing potential, respectively. The absorbing potential should be sufficiently 
strong to absorb all the outgoing waves, and should be smooth enough to avoid the occur- 
rence of the significant reflected waves. These conditions could be fulfilled if we employ a 
large enough value for Ar. The large Ar, however , implies the necessity of treating a large 
spatial area, which in turn results in the increase of the computational task. In practice, the 
following condition derived in the WKB approximation, gives a criterion for a good absorber 
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20—-= < W^abs < ■77:^r-^ . 8 

ArV8m 10 ^ 

The left inequality of Eq. (^) originates from the condition that the absorption is strong 
enough to suppress any reflection at r = rc -|- Ar, while the right inequality originates from 
the condition that the reflection at r = r^ is sufficiently small. Lower the bombarding energy 
E becomes, the wider Ar is required to fulfill the condition of Eq. (||). This is due to the 
increase of the de-Broglie wave length at low energies. 

In order to demonstrate the quality of the ABC calculation, we consider, as an example, 
the elastic scattering of ^^O + ^^C at i^iab. = 139.2 MeV with the optical potential of 



Ref. [|T3l. The potential V{r) is spherical and composed of the real and imaginary parts of 
the nuclear potential and of the Coulomb potential whose short range part is regularized 
with a quadratic function (the Coulomb potential of a uniformly charged sphere). 

Although one can easily obtain the numerical solution for the potential V"(r), we divide 
the potential into two parts, Vo(r) and AV"(r), in order to test the ABC. Here, we adopt the 
real part of the nuclear potential as AV{r), and the rest as Vo(r). We solve Eq. (^) in the 
partial wave expansion, Aip'^^^r) = I]/m(2//(^)/'")^m(f")- 

We fix the parameters of the absorbing potential, Eq. (|^), as = 20 fm and Ar = 10 fm, 
but vary VFabs- The vanishing boundary condition, yi{r) = 0, is imposed at r = rc + Ar = 30 
fm. The discrete variable representation, which will be explained in detail in Section III, is 
utilized with the uniform grid spacing of 0.2 fm. 

In Fig. (a), we show absolute values of the 5- matrices as a function of I calculated with 
different values of W^bs ■ The closed circles represent the absolute values of the S matrices 
which are obtained by a standard numerical procedure to integrate the radial Schrodinger 
equation up to r = 30 fm. One can see that the ABC method works quite well in describing 
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FIG. 1: (a) The absolute values of the S'-matrices as a function of the relative angular momentum 
/ for the elastic ^^O + ^^C scattering at -ELab. = 139.2 MeV. (b) The real part of the radial wave 
function of the scattering partial wave / = 30. 
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FIG. 2: The elastic scattering S matrix elements Si (a) and the differential cross section (b) for 
the same reactions as shown in Fig. |l|. 

the potential scattering. It should be noted that only the absorbing potential with IVabs = 
200 MeV satisfies the condition of Eq. (|). 

In Fig. [l| (b), we show the real part of yi{r) for the partial wave of / = 30, which 
corresponds to the grazing angular momentum. The three lines correspond to the different 
choices of the strength IVabs- The damping behavior of the three wave functions in the region 
20 < r < 30 fm are different according to the strength Wabs- Despite of small deviation 
in the wave functions at r < 20 fm, it does not affect the absolute values of the elastic 
scattering S matrix elements as we saw above. This implies that the reaction cross section 
is not so sensitive to the choice of the absorbing potential. 
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However, when evaluating the differential cross sections, one has to be careful on the 
choice of the absorbing potential. Figure ^ displays the Argand diagram of the elastic 
scattering S matrix elements (a), and the differential cross sections (6). The elastic scattering 
cross sections in the Rutherford ratio vary by a few orders of magnitude as a function of 
the scattering angle. The subtle difference in the scattering matrix results in the substantial 
deviation in the differential cross sections. The calculation with pyabs=200 MeV, which 
fulfills the condition of Eq. (|]), is reasonably accurate while the other choices of W^hs, which 
do not fulfill the condition, give less accurate results at large angles. 



III. ABSORBING BOUNDARY CONDITION METHOD: A THREE-BODY RE- 
ACTION 



A. Basic equations 

In this section we explain how to apply the ABC to three-body breakup reactions. We 
consider the reaction of a projectile (P) on a target (T), in which the projectile is composed 
of a core nucleus (C) and a neutron (n). We denote masses of the core, the neutron, and 
the target as m^, m„, and m-p, respectively. The charge numbers of the core and target are 
Zc and Zt, respectively. We employ a Jacobi coordinate system in which the n-C relative 
coordinate is denoted as r and the P-T relative coordinate as R. We assume that the total 
Hamiltonian is given by the following form, 

H = -^V^ - —Vl + V;c(r„c) + KT(r„T) + Vcr(rcT), (9) 

where the reduced masses, /i and m, of the relative P-T and n-C motions, respectively, are 
given by 

^ = 7 r— , m = • (10) 

(m„ + mc) + mT rrin + mc 

The vectors 

r„c = r , r„T = R H \ r , tct = R \ r, (11) 

rUn + mc mn + mc 

represent the n-C, n-T, and C-T separations, respectively. Vnc, VnT, Vct are the effective 
potentials between the corresponding particles. Vct includes the Coulomb potential as well 
as the nuclear potential. VnT and Vct may be complex. Their imaginary parts represent 
excitations of the target and/or core nucleus. Vnc is assumed to be real. 



We introduce an auxiliary distorting potential Upt(R) for the P-T relative motion. The 
Upt(R) niay be chosen arbitrarily as long as it contains the long-range Coulomb part of 
the P-T interaction. The distorted wave function $k*^^'' (R-) is a scattering solution of the 
following equation with the outgoing boundary condition, 
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-|;V^ + C/pt(R) 



$S(-)(R) = £;$£W(R), = (12) 

where K is the incident wave number vector. In the partial wave expansion, we express 
$^^+Hr) as 

$g(+)(R) = 47rEi^e^(-+^-):^^i^y;^(K)y^M(R). (13) 

For the n-C motion described with the Hamiltonian hnc — + Kic(r), we denote 

the solutions as 0/3(r) , where /3 specifies the quantum numbers n^l^mfj for the bound 
orbitals and kplpmp for the continuum states, np is the nodal quantum number and kp is 
the magnitude of the wave number. Ip and mp are the orbital angular momentum quantum 
numbers as usual. In the polar coordinate system, 0/3 (r) is rewritten as 

0,(r) = ^y,^„^(f), (14) 

where up — Ungig for the bound orbitals and Ujs — Ukgig for the continuum states . The 
bound wave function is normalized as usual, and the continuum wave function is normalized 
by the following asymptotic form, 

fJ-k/shir) ~ y^si" (^h'^ ~ '^'z') ' r ^ oo, (15) 

where Si^ is the phase shift for the n-C scattering. 

We express the total wave function ^Koq(R, r) as the sum of the distorted wave in the 
incident channel a and a scattered wave as 

<l(R,r) = <+nR)</'a(r) + *s(R,r), (16) 

where Kq is the incident wave number vector of the P-T relative motion. 
The scattered wave funtion ^s(R, r), satisfies the following equation 

{E - H) *s(R, r) = Al^(R, r) (R) ^.(r), (17) 

with 

A1/(R, r) = V^tM + Vct{vct) - Upt{K) , (18) 



where E is the total energy which is equal to the sum of the energies of the incident P-T 
relative motion h'^'Kl/2fi and the energy of the n-C relative motion in the incident channel, 
The right hand side, AV^(R, r)$^[|^^(R)0a(r), should be localized in space. In other 
words, this function vanishes if either or r is larger than a certain critical radius, Rc for 
R and for r [^. 

Since the scattered wave function, \I's(R-, i"), includes only the outgoing waves asymptot- 
ically, it should be legitimate to simulate the outgoing boundary condition by introducing 
the absorbing potentials for both R and r coordinates. Thus, instead of solving Eq. {\[7\] 
with the outgoing boundary condition, we propose to employ the following equation with 
the vanishing boundary condition. 

[E + lencir) + iepT{R) - H] ^s(R, r) = Al^(R, r) ^g^+^R) 0a(r), (19) 

where enc{.f) and epT^R) are the absorbing potentials for the relative motions of n-C and 
P-T, respectively. They are placed in the spatial region R > Rc for epriR), and r > 
rc for enc{^)- In practice, we employ the linear absorbing potentials of Eq. (|^) for both 
functions. The vanishing boundary condition is imposed for \E's(R, r) at -R = Rc + Ai? and 
r = r(7 + Ar. In this way, the three-body scattering problem is converted into the three-body 
Schrodinger-like equation with the vanishing boundary condition. No explicit consideration 
on the boundary conditions for three-body continuum states is required in this approach. 

As for the parameters of the absorbing potentials, careful considerations are required. For 
iepT^R), we choose the potential parameters according to the condition of Eq. (H) , where 
the energy E may be evaluated at the relative P-T energy in the incident channel. However, 
for i€nc{f^), the energy of n-C motion after reaction is not unique but spreads over a certain 
energy region. A reasonable choice is then to set the parameters of ienci^) optimum for 
the main component of the n-C motion. It should, however, be noted that the complete 
absorption of the breakup waves is difficult for small n-C relative energies because of their 
long wavelengths. 

The three-body wave function obtained in the above procedure is meaningful in the inner 
spatial region, r < rc and R < Rc- The wave function only in this interacting region is 
required to calculate the relevant T-matrices. For the breakup processes in which the final 
n-C motion is specified by the relative wave number k and the final P-T motion specified 
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by K, the relevant T-matrices are given by 

T,(K,k) = {^^^-\R)(^i-\r)\VnT + VcT-UpT\^^^liR,r)) , (20) 

where (f)^\r) is the n-C wave function with the incident plane wave specified by the wave 
number vector k and with the incoming boundary condition. This is related to ipkimi^) of 
Eq. dll by 

<pVij) = E ^^^'e-'^'0,,^(r)l^:;(k) , (21) 

Im 

where 6i is the phase shift for the n-C scattering. For convenience, we introduce an alter- 
native definition for the T-matrices, Tq,/q,(K), in which a' specifies the final n — C state, 
a' = nlm for the bound orbitals and a' = klm for the continuum states. This T-matrix is 
given by 

T„,,(K) = (<l>g(-)(R)</.,,(r)|KT + ^CT-f/pT|^£jR,r)). (22) 

For the reactions in which the final state is the three-body continuum state specified by 
a' = klm, the T-matrices of Eqs. (pO]) and (^) are related by 

T„(K, k) = ^(-0'e^^'y;^(k)Tw^,«(K). (23) 

Im 

Various cross sections can be expressed in terms of these T-matrices. For example, the 
angle-integrated elastic breakup cross section is given by 

\27rh / Ao 
B. Radial coupled channel equation 



In practice, Eq. (|19| ) is solved in the partial wave expansion. The incident wave function 
•^Ko (R')0a(r) and the total wave function \l/^j^(R, r) are expressed as 

$g(R)0«(r) = \/4^(2^ + l)(i^O/,m,| Jm,)z^e^(-^+'5^) 



JL 



X ^^^f|^^^[r.(R)rjr)],.„, (25) 



and 



^Kol(R,r) = \/47r(2L + l)(L0/«m,|Jm«)^^e*('^^+'^) 



JLL'V 



X , (26) 
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respectively . For the scattered wave function \1/5''^(R, r), we make the same expansion as 
Eq. (p6D , replacing the radial wave function vii^ Liii{R,r) with aj^i^ ^,i,{R,r). These two 
radial wave functions, y[i^ j^iii{R,r) and aj^i^ j;^,i,{R,r), are related to each other by 

1 



yLi^,L'i'iR^r) = — fL^Ko, R)un^iSr)hL'Ki' + (^Li^,L'i'iR^r). 



(27) 



The radial coupled-channel equation for the scattered wave function {R, r) is writ- 



ten as 



E + itncir) + iepT{R) 

j n^L'{L' + i) d hH'(i' + i) 



+ 



- ^L'v,L"i"i.R^'r)aii^^L"U'i.R^r) 



+ Vncir 



a 



J 

Lin 'L^ I 



{R,r) 



L"l" 

= [vi.,^uSR.r) - hv^i^vUpTiR)} ^JdK.^RW^ar) , 
where the radial coupling potential Vlii, j^„i„{R,r) is defined by 

Vln>u'V'{R.r) = I f/Rrfr[FL'(R)>^z'(i-)]}M(Kr + VbT)[lL.(R)ll.(f)] 



JM 



(28) 



(29) 



The T-matrix elements are expressed in terms of the radial wave functions . For a 
transition concerning the n-C relative motion from the bound orbital a = riala^a in the 
incident channel to the continuum state a' = klm in the final channel, the relevant T-matrix 
element is expressed as 

T,,,(K) = {^^^-\n)Mv)\VnT + VcT - UpT\^^^l^{n,v)) 

= ^ (47r)^''^(2L + l)l/2^i-i'gi'T£+ia^,+i5i,+j<5^, 
JLV 

x{LOlama\Jma){L'ma - mlm\Jma)YL'm^-rn{^)li,L,^aL ) (30) 

where the radial integral I'^ij^i ^.^ is given by 

li'L',aL = E / dRdr-hiK, R)u^.{r) {vI,i^l"1"{R. r) - 6l'l"Su»Upt{R)} yii^,L"i"{R, r). 

(31) 

For the elastic channel in which the final state a' is identical to the initial one, a' = a = 
rtala^ai the two-body T-matrix for Upt should be added to the above expression of Eq. (pO|). 

For later convenience, we present an expression for the angle-integrated elastic breakup 
cross section in terms of this radial integral. The elastic breakup cross section can be 



11 



decomposed into the sum over the total angular momentum J and the relative angular 
momentum / of the n-C motion, as well as the integral over the relative wave number k, 



poo 

^breakup ^ ^ / dkaU{k), 



with 



aU{k) = levr 



2 K 



L 



(32) 



(33) 



where K is the P-T relative wave number which implicitly depends on the n-C wave number 
k through the energy conservation relation. 



C. Computational details 

To solve the radial coupled-channel equation (0), we discretize the radial variable with a 
uniform grid. To achieve high accuracy with a moderate number of grid points, we employ 
the so-called discrete variable representation (DVR) [1^. Although the DVR is a kind of 
basis expansion methods, the resultant equation is similar to that in the finite difference 
approximation. 

We employ labels i and j for the grid points of the radial coordinates R and r, respectively. 
Denoting the grid spacings for the coordinates R and r as H and h, respectively, the grid 
points Ri and rj are given by 

Ri = tH ,{t = l---NR) rj=jh ,{j = 1---Nr) , (34) 

where Rc + AR = Nji ■ H and + Ar = ■ h. The wave functions are then discretized as 

In the DVR, the matrix elements of the potentials become diagonal, 

\^Ll,L'v]i3,i'j' = ^Ll,L'v{Rh'^j)^ii'^j3' ■ (36) 

The kinetic energy operator for the coordinate R is expressed as the following matrix form, 

2/idi?2 



5 ^ D« 



(37) 
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where Df^, is defined by 

r 2^ _ ^ for = t') 

(38) 



for [i ^ i!) 



Tlie kinetic energy operator for the radial coordinate r is expressed in a similar manner. 

After these procedures, the radial coupled-channel equation, Eq. (PBD, is expressed as the 
following matrix form, 

XI [^L'i',L"/"]ii,i'i'[«U„,L"Z"]i'i' = [^L«a,L'«']ii > (39) 
L"«"i'j' 

where the coefficient matrix [A'l'i\L"i"]ij,i'j' is defined by 



[^L'l',L"l"]ij,i'j' ~ E5L'L"5l'l"5ii'6 



(Ti^L'{L' + l) \ 1 

+ ( 1^^2 iepT{Ri)\ 5ii>\5jji5L'L"Svi" 



rj 



-Vuv,u'i"iRi^ rjW'Sjf (40) 
and the source term [bLi^^L'i']ij is given by 

[bLl^,L'vh = {Vli^^LiA^i^ ''j) - hL'Sll^UpT{R^)} ^/l(^0, RiW^l^iVj) . (41) 

The dimension N^i^ of this linear algebraic equation is given by the product of the number 
of angular momentum channels , which is denoted by Nj, the numbers of grid points A^^^; and 
Nr , i.e., Ndim = NjNuNr. As will be seen in a practical example given in the next section, 
this dimension becomes as large as A"dim = 10^. The coefficient matrix [Aj^ii, L/'i'/]ij^i'j' is 
sparse, namely, the number of non-zero elements of this matrix is rather small. 

There are many efficient algorithms which are proposed to solve such a linear algebraic 
equation with a sparse coefficient matrix. We will employ the Bi-conjugate gradient method 
]T3| , which is one of the most well-known iterative method to solve the linear algebraic equa- 
tion with a complex non-hermitian coefficient matrix. In order to improve the convergence, 
the pre-conditioning [|l6l is carried out. In our problem, the diagonal part of the matrix A 



includes huge elements because of the centrifugal barrier. In order to balance the matrix, 
we modify the equation as 

^ A^ • VBd = -^h , (42) 
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where D is the diagonal matrix whose elements are equal to the diagonal matrix elements 
of A. The convergence of the Bi-conjugate gradient method is significantly accelerated for 
the matrix -^A^^ compared with the original matrix A. 

After solving the equation, we calculate the T-matrices with the similar discretization. 



IV. DEUTERON REACTIONS: COMPARISON WITH THE CDCC METHOD 



We here apply the ABC method to the rf+^^Ni reaction at the incident deuteron energy 



Ed = 80 MeV. Since this reaction was studied in detail with the CDCC method fl^ , we can 
assess validity of the ABC method by comparing our results with those of Ref. |1T7| . 

We employ the same Hamiltonian as used in Ref. Namely, the global optical potential 
for nucleon-nucleus elastic scattering [0 is adopted as the interactions between the proton 
(neutron) in the deuteron and ^'^Ni. The potentials at half of the incident deuteron energy 
is adopted. No spin-orbit force is taken into account. The potential parameters are: Vr = 
44.921 MeV, tr = 1.17 fm, aR = 0.75 fm, Wy = 6.100 MeV, Wsf = 2.214 MeV, rj = 1.32 
fm, ai = 0.534 fm for proton and Vr = 42.627 MeV, = 1.17 fm, qr = 0.75 fm, Wy = 7.240 
MeV, Wsf = 2.586 MeV, r/ = 1.26 fm, aj = 0.534 fm for neutron. The Coulomb potential 
of uniformly charged sphere is assumed between the centers of mass of the deuteron and 
^^Ni. The radius of the sphere is set as 4.92 fm. 

The potential for the p-n relative motion is taken as a Gaussian form, Vnp{r) = 
Voexp[— r^/ag] with Vq = —72.15 MeV and = 1.484 fm. The ground state is assumed to 
be a pure s-wave. This potential gives the deuteron binding energy of 2.22 MeV. 

We treat the radial regions of R and r up to 50 fm. The radial scattered wave function 
(y[i^ Liii {R, r) is set equal to zero if either i? or r is larger than 50 fm. The absorbing 
potentials iepxiR) ^e„(7(r) are active at R,r > 25 fm. They are parameterized as 

tepT{R) = 'WpT^^, ^enc{r) = ^"^nJ-^ , (43) 

with WpT = 50 MeV, R^ = 25 fm, Ai? = 25 fm, Wnc = 20 MeV, = 25 fm, and Ar = 25 
fm. This implies that the obtained solution will be reliable in the spatial region oi R,r < 25 
fm. As for the auxiliary distorting potential Upt(R.), we take the folding potential in which 
the p-T and n-T potentials are folded with the deuteron ground state density. 

We carry out the calculations for total angular momentum J up to 100, including / = 



14 



(a) L' = 17:l' = 



(b) L' = 15:l' = 2 




FIG. 3: (a) The real part of the radial wave function y^^^ I'l'i^^ ^) ^'^^ ~ ^"^^ ^' ~ ^"^^ /' = 
in the d-^^Ni scattering at Ed = 80 MeV. (b) The same as (a) but for J = 17, L' = 15, and = 2. 

and 2 for the p-n relative angular momentum. The radial discretization is made with a 
spacing oi H = 0.2 fm for i? and /i = 0.5 fm for r. The matrix size is then given by the 
product of the number of angular momentum channel Nj = 4, the number of radial grid 
points Nji = 250 and Nr = 100. This amounts to the total number of points N^i^ = 10^. 
As mentioned in the previous section, we solve this linear algebraic problem with the sparse 
complex matrix using the Bi-conjugate-gradient method. 

As for the grid spacing of the p-n coordinate, we found the spacing of 0.5 fm gives the 
accuracy of the deuteron binding energy better than 0.0001 MeV. The accurate binding 
energy by the n-p potential employed is 2.2177 MeV, while the DVR calculations with the 
grid spacing of 0.1 fm, 0.2 fm, and 0.5 fm give the same value of 2.2177 MeV. 

In the present approach, we can explicitly obtain the wave function. In Fig. ^ (a), the 
real part of the radial wave function y[i^ j^,i,{R,r) is displayed ioi J = L = L' = 17, and 
I' = la = 0. This includes the elastic scattering wave as well as the breakup waves. In the 
small r region the wave function is dominated by the elastic wave . It oscillates with the 
frequency close to Kq in the coordinate of R. Its amplitude decreases as R increases in the 
region R > 25 fm, because of the absorption by —iepT{R)- The wave function at large r 
shows breakup components of the deuteron into the p + n continuum state. The amplitude 
of the wave function also decreases at large r due to the absorption by —ienc{i"). Fig. |^ 
(b) shows the real part of viicVvi^^''^) — 1''') ^' — ^5, and /' = 2. For this angular 

momentum channel, only the breakup component is seen. 

We next show decomposition of the breakup cross section cr/^(fc), which was defined in 
Eq. (|33D . Following Ref. Jl^, we show, in Fig. ||, the square moduli of the 5*- matrix, which 
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FIG. 4: (a) The squared moduli of the deuteron breakup 5-matrix elements |5^^j;^Qp. (b) 
|5^^j_2;=2p- (*^) \S'i=ji=2\'^ ■ (d) \S'[^jj^2i=2\'^ ■ each panel, the dotted line and histogram 
denote the results of the ABC and CDCC methods, respectively. 

are related to cr'lAk) by crl,(k) = ^^^^^P^\Sii\k)\'^ for the present case (la = 0). The deuteron 
s-wave breakup, S'{{^j^^q(A;), are shown in the panel (a) for J = 17. Three kinds of rf-wave 
breakup S"- matrix elements are shown in the panels (b) S^j^2j_2i=2{k)> (c) S^Llji=2{k) ■> ^'^^ 
(d) S'flj^2i=2{^)i respectively. In each panel the dotted line and histogram represent the 
calculated results of the ABC and CDCC methods, respectively. The global behavior of the 
smooth dotted lines agrees well with that of the histogram . We examined the sensitivity 
of the breakup cross sections against the strength of the absorbing boundary potential. We 
confirmed that the calculated results are not sensitive to these parameters, as long as the 
condition of Eq. (|) is satisfied. 

We show in Table | the total reaction cross section aR and the elastic breakup cross 
sections asu ■ The elastic breakup cross sections are decomposed according to the p-n 
relative angular momentum /. While the total reaction cross section and the s-wave 
breakup cross section 0"^^°'' calculated with the ABC method agree to those with the CDCC 
method, the rf-wave breakup cross section 0"^^^'' of the ABC method is about 10 % smaller 
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TABLE I: Reaction and breakup cross sections of d- Ni reaction at Ed = 80 MeV. The results of 
CDCC are taken from Ref. 0. 



(/=0) (1=2) 
'^BU '^BU ^BU (TR 

ABC 38.28 mb 77.21 mb 115.49 mb 1572.9 mb 
CDCC 37.77 mb 88.59 mb 126.36 mb 1571.4 mb 
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FIG. 5: The S'-matrix (a) and differential cross section (b) of d-^^Ni elastic scattering at Ed = 
80 MeV. Results of the ABC are compared with those of the CDCC [||. Results using a folding 
potential are also shown in (b). 

(1=2) 

than that of the CDCC method. This difference in a^^ is attributed to the difference seen 
in Fig. |(b). 

We next consider the angular distribution of the elastic scattering. Fig. ^ (a) shows 
the elastic scattering S'-matrix elements Sj. The closed circles with dotted line denotes Sj 
calculated with the ABC methods. The open squares represent the corresponding CDCC 
results for J = 1, 8, 12, 15, 16, 18, and 20. Fig. || (b) shows the elastic scattering differential 
cross section in Rutherford ratio. The differential cross section calculated with the folding 
potential is also shown. The difference between the ABC (CDCC) calculations and the 
folding potential calculation shows the effect of the breakup process on the elastic scattering. 
The results clearly show that the ABC method gives almost the same results as the CDCC 
results. Therefore, the ABC method can describe the influence of the breakup processes on 
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the elastic scattering angular distribution. 



V. SUMMARY 

The purpose of the present article is to show the usefulness of the ABC method to describe 
nuclear breakup reactions. The basic trick in the method is to place an absorbing potential 
in the outer spatial region where complex three-body reactions do not take place. The 
vanishing boundary condition is imposed outside the absorbing potential. If the absorbing 
potential works properly, there exist only the outgoing waves at the edge of the interacting 
region where the absorbing potential is absent. We can thus obtain scattering solution by 
solving the Schrodinger-like equation with a source term, without imposing the scattering 
boundary condition. 

In practical calculations, the equation is solved in the partial wave expansion, and the 
radial equation is discretized employing a uniform grid. The problem is then converted into 
the linear algebraic equation of large dimension with a complex, sparse, coefficient matrix. 
We can efficiently solve the equation by employing the Bi-conjugate-gradient method with 
a pre-conditioning procedure. 

To show the feasibility and the accuracy of the ABC method, we have applied the method 
to the (i-^^Ni scattering as an example. For this scattering, detailed analyses with the CDCC 
method are available. We have found an excellent agreement between the results by the 
two methods, the ABC and the CDCC methods, for both breakup and elastic scattering. 
Prom this observation, we are confident that the physical contents of these two methods are 
the same, although the procedures in solving the three-body scattering problem look very 
different. 

Let us compare the advantages and the limitations between the two methods. One of the 
superior features of the ABC method is that it is easy to assess convergence of the calculated 
results. If we fix the radial grid points and the number of partial waves for the n-C motion, 
the quality of the absorbing potential determines accuracy of the results. One can easily 
ascertain the reliability by examining how sensitive the results are to changes of the absorbing 
potential. The convergence check is more difficult in the CDCC method in which one must 
ffist prepare the discretized continuum wave functions. One must examine the delicate 
convergence as to the discretization of the continuum; the number of the momentum bins 
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and the maximum momentum to be incorporated. 

The principal drawback of the ABC method is its heavy computational cost. One must 
solve a linear algebraic equation of large dimension. In the present example of deuteron 
reaction, one must solve the equation of 10^ dimension for each angular momentum. For 
this problem, the computational time for each angular momentum is typically 20 minutes in 
a single alpha-21264 processor. This time multiplied by the number of angular momentum 
states gives the total computational time for the reaction at a fixed incident energy. Of 
course, the parallel computation would reduce the computational wall-time, by paralleling 
different angular momentum calculations. 

In general, as the incident energy gets higher, the calculation becomes more demanding. 
This is because the finer grid spacing is required as the incident energy becomes higher. 
Also the number of angular momentum states that take part in the reaction increases. For 
reactions with heavy target nuclei, the Coulomb breakup reactions are significant, in partic- 
ular when the projectile is bound extremely weak, such as the halo nuclei. The description 
of such Coulomb breakup reactions will be a difficult problem for the ABC method, because 
the coupling potential is so long-ranged that one should employ a large critical radius Vc to 
obtain reliable results. We must place the absorbing potential at a large radial distance, and 
must employ a large number of grid points. These difficulties are, however, not inherent in 
the ABC method, but are suffered by the other methods as well. 

Establishing the usefulness of the ABC approach for breakup reactions of a weakly bound 
projectile, we are now planning to apply the method to reactions of unstable nuclei. A 
preliminary result on the nuclear breakup reactions of a single neutron halo nucleus, ^^Be, 



will be published in Ref. pO[. In that calculation, the number of the matrix size is about 
2.4 X 10^, not so much different from that in the deuteron reaction. The ABC method 
should, at least, be useful for breakup reactions of a single nucleon halo nucleus induced 
by the nuclear force. Applicability of the ABC method depends on the availability of the 
computational facilities. We expect that the ABC method will be more popular and useful 
as the parallel computation becomes more prevalent. 

The ABC method is also useful in any circumstances where the scattering boundary 
condition comes in. Some applications have been done and are in progress concerning 
resonances of a few-body system and response in the continuum, especially, of deformed 



systems [12, 21, 22 
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